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A THEORY FOR mEDICTIRG COMPOSITE LAMINATE 
WARPAGE RESULTING FROM FABRICATION 

By C . C. Chamis* 

NASA-Lewis Research Center 
C le ve land , Ohio 

ABSTRACT 

Linear laminate theory is used in conjionction with the moment- curvature 
relationship to derive eqviations for predicting end deflections due to war page 
without solving the coupled fourth-order partial differential equations of the 
plate. Conposite micro- and macromechanics are \ised in conj\mction with laminate 
theory to assess the contribution of factors such as ply misorientation, fiber 
migration, and fiber and/or void volume ratio nonuniformity on the laminate war- 
page. Using these equations, it was found that a l'^ error in the orientation 
angle of one ply was sufficient to produce warpage end deflection eqml to two 
laminate thicknesses in a 10 inch by 10 inch laminate made from 8 ply Mod-l/epoxy. 
Using a sensitivity analysis on the governing parameters, it was found that a 3 ° 
fiber migration or a void volume ratio of three percent in some plies is sufficient 
to produce laminate warpage corner deflection equal to several laminate thicknesses. 
Tabular and graphical data €o:e presented which can be used to identify possible 
errors contributing to ’ 'jninate warpage and/or to obtain an a priori assessment 
when unavoidable errors during fabrication are anticipated. 
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INTRODUCTION 


Flat composite laminates have been observed to warp upon removal from their 
fabrication mold. This warpage is the result of nonsyrametries and residual 
stresses that are present in the laminate. The warpage can be of a magnitude 
sufficient to render the laminate useless for its intended purpose. Thus, residxial 
stresses are always present in angleplied laminates which are cured at elevated 
temperatures and then cooled down to room temperature. Bending nonsymmetries 
result from ply misorientations and fiber nonuniformities which occur inadvertently 
during the fabrication of the laminate and produce coupling modes in the laminate. 

When the laminate is subjected to either thermal or mechanical loads, these 
coTipling modes produce laminate warpage. The various co\ipling modes resulting 
fi*om combinations of nonsy uimeti' i e s are discussed in reference 1. Exact deteinnina 
tlon of the surface of the warped l am i n ate requires solution of coupled fourth 
order partial differential equations. In the literature, neither solutions to the 
coupled partial differential equations or approximate equations are available for 
predicting laminate warpege as a result of bending nonsymmetries and residual 
stresses. Therefore, the assessment of the effect of various factors that contri- 
bute to the warpage and the establishment of fabrication control procediires has not 
been possible in the past. 

Equations, even of an approximate nature, for predicting laminate warpage 
could be of considerable practical value to both fabricators and designers. For 
exan 5 )le, the fabricator could use the equations to obtain an a priori assessment 
of the warpege resulting from factors which may be difficult to control accurately 
during laminate fabrication. The designer, on the other hand, could use the results 
to specify tolerances which woiild minimize warpage-producing nonsymmetries, or he 
could use the results to design the laminate with warpage-compensating nonsymmetries. 
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It Is the objective of this report to rhow how linear laminate theory can 
be used to derive convenient equations without solving the coupled fourth-order 
partial differential eq^oations as would normally be required. The ase of these 
equations for predicting warpage is demonstrated. 

WARIAGE GE(»ffiTRY AND ORIGIN 
Warpage Geometry 

The laminate warpage geometry of interest in this investigation is depicted 
in figu"*. This type of warpage is typical of \ laminates which have 
the 0° ply direction parallel to the ::-axis. As can be seen in figure 1, the laminate 
is fixed along its x-edge (AB). Note that in this position, the laminate exhibits 
two primary modes of warpage. These are: (l) curvature (bending) along the y- 

direction and (2) a twist about an axis normal to and bisecting the plate x-edge. 

There is also a curvature along the x-direction. However, for the laminate con- 
figurations under consideration, the curvature along the x-direction is relatively 
small compared to that of the y-direction. The corner deflections C and D provide 
a quantitative measure of the amount of warpage. The twisting increases the bending 
corner deflection at C while it decreases the corner deflection at D, In an actual 
warped Itminate, the corner deflections C and D can be readily measured by orienting 
the laminate as is shown in figure 1. 

Warpage Origin 

Composite laminates are, in essence, nonuniform materials since they consist 
of several plies with different orientations through the thickness. Becaiise of this 
nonuniformity and the lay-up procedure used diuring fabrication, some bending- stretching 
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coupling in the laminate is the general rule rather than the exception. Bending- 
stretching coupling is present when a purely extensional deformation in the laminate 
produces, simultaneously, bending and stretching. Conversely, a purely bending or 
twisting deformation in the l€uninate can produce, simultaneously, an extensional 
deformation. Bending- twisting coupling is also possible in angleplied laminates. 

This is the case even if the angleplied laminate is balanced and symmetric. 

Both ply misorientations and fiber nonuniformity will produce bending- stretching 
coupling which resvilts in warpage of a flat laminate in the presence of residual 
stress. When these conditions are present in laminates, they will 

produce warpage having the geometry depicted in figure 1. 

THEORETICAL BACKGROUND AND GOVERNING EQUATIONS 

In this section, the theoretical background leading to the governing equations 
for predicting laminate warpage corner deflections in the presence of residual 
stress are described. The equations presented for approximate deflection prediction 
are of convenient form and obtainable without solving the coupled fourth-order 
partial differential equations. A list of symbols used in the equations is con- 
tained in the Appendix of this report. Warpage corner deflections can be predicted 
more accurately by anisotropic thin-plate bending theory whi.ch accounts for all 
possible cotipllng responies present in such plates (Appendix, reference l). Warpage 
corner deflections can also be predicted accurately by special finite elements 
accounting for the same coupling responses. However, neither anisotropic thin- 
plate bending theory nor finite element analysis which accounts for all coupling 
responses have been used to predict laminate warpage corner deflections in the 
presence of residual stress. 
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Theoretical Background 

The derivation of the equations for predicting laminate warpage corner 
deflections is hesed on the force-deformation-temperature relationships derivable 
from linear thin- laminate bending theory (refs. 2 and 3). The two important assump- 
tions in linear thin-laminate bending theory are: (l) the planform dimensions of 

the laminate are very large relative to its thickness and (2) the maximum deflection 
of the laminate is of the same order of magnitude as the laminate thickness. The 
second assumption is usually interpreted to mean accurate predictions for maximum 
deflections equal to plate thickness and good approximation for maximum deflections 
of 10-20 times the plate thickness. The 10-20 times facuor is based on the maximum- 
deflection to plate-edge— dimensions ratio. If this ratio is such that the sine and 
tangent of man angles can be approximated by the angle itself (in radians) then 
thin-laminate bending theory would predict good approximations to maximum 
deflections . 

To illustrate this point with an example, consider a (0g,^5,+^5,02) laminate 
with m»Tr-iTnmit tip diflection of 1.20 inches, plate edges of 10 inches, and plate 
thickness of .06 inches. The tangent of the angle subtended by the 1.20 inch 
dimension on a 10 inch base is 0.12. The corresponding angle in radians is 0.119 
^fhich is equal tr the tangent when rounded off to two figures. Based on the afore- 
mentioned criteria, linear thin laminate bending theory would be applicable to these 
types of laminates even though they undergo corner deflections 20 times their laminate 
thickness. 

The laminate force deformation relationships derivable from thin-laminate 
bending theory and of interest in this investigation are given in references k or 
5. In matrix form and assuming forces due to residual stresses only, these equations 
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become: 
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(1) 


The notation in equation (l) is as follows: denotes reference plane 

strains denotes the plate local curvatures; A represents a(3 x 3) array 

cx 

of axial (membrane or stretching) stiffness coefficients; represents a (3 x 3) 
array of stretching-bending coupling stiffhess coefficients; represents a 

(3x3) array of bending (flexural) rigidities coefficients; denotes themml 

forces due to residual stresses and are the corresponding thermal moments. 

Note the subscript x in eqx^tion (l) indicates that these relationships are referred 
to the laminates’ structural axes (x,y,z, fig. l). Note also the superscript -1 
denotes the inverse of the matrix. 

The equations used to generate the elements in the arrays A , C .and 

^ cx’ cx ’ 

^cx given in reference 4. The form of these eq;mitions, neglecting interply 

contribution, is: 

NZ 

= ^U(^zi+1 - (2) 

i=l 

[Ccxi = £ (3) 


'°cx'=C 

i=l 




(4) 
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The corresponding equations for the thermal forces and moments are: 
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The arraysSl^^i}, equations (2)-{6) are: 
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(5) 

( 6 ) 


( 7 ) 


( 8 ) 


(9) 


The notation in equations (2)-(9> is as follows. R, is the number of pUes 
in the laminate; is the temperature of the i-th ply or the difference be- 
tween cure and room temperatures in the present case; the difference - 

2^. locates the i-th ply relative to the reference plane (stacking sequence); 
o^lll is the ply thermal coefficient of expansion along the fiber direction and 

normal to it; is the angle locating the fiber direction in a ply relative 

to the laminate structural axes (ply orientation angle); is the ply modulus 

of elasticity along the fiber direction and E^^f normal to it; is the major 
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Poisson’s ratio and 5^11 inplane ply shear mod\ilus. 

It is well known that the ply thermal coefficients of expansion, moduU and Pbisson’s 
ratios are related to constituent material properties using composite micromechanics 
when the types of constituents, the fiber volume ratio, and the void volume ratio 
are known. See reference 4. 

The variables of interest in the present discussion are the laminate local 
curvatures {X.c% ), defined in equation (l) which contribute to laminate warpage. 

It is seen hy inspection from eqtiations (l)-(9) that the laminate curvatures depend 
on: the nunfiber of plies in the laminate (Nj ), the ply t«nperature difference 

(A.T), the ply stacking sequence ), the ply orientation angle 

( Oji, ), and through micromechanics on the constituent material properties (E, G, 

V, ando<), the fiber volume ratio, and the void volume ratio. 

Equations (l)-(9) and the micromechanics equations for predicting ply proper- 
ties from constituent properties have been programmed in the computer code described 
in references 4 and 5. This computer code is used herein to compute the laminate 

local curvature JC explicitly as a function of the factors identified in the previous 
OX 

' iph. Therefore, it becomes strai^t forward to compute perturbations of the 
fact'is contributing to warpage about the condition of bending syirmetry. How the 
local curvat\ires will be used to predict laminate warpage corner deflections will 

be described in the second part of this section. 

The three important points to be noted from the previous discussion are. 

1. Laminate waipage, of interest in the present discussion, is within the 
realm of linear thin- laminate bending theory even when the laminates 
undergo large corner deflections compared to their thickness. 
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2. The factors contributing to laminate wax'page are readily identified via 
the laminate force-defonnation-teniperature relations. 

3. A cmtputer program is available that csin be used to assess the relative 
significance of the various factors contributing to laminate warpage. 

Governing Equations 

The governing equations for predicting laminate warpage corner deflections 
are readily derived when the local curvatures are Icncwn and advantage is taken of 
the following observation. The laminate local curvatures (l)) resulting 

from uniform residual stress ( ^T constant in eqs, (2)-(6); are constant throu^out 
the laminate. This observation is violated in the vicinity of the free edges within 
a distance apprCTcimately equal to the laminate thickness which is Insignificant 
coaqwred to laminate planform dimensions. Since the radius of curvature is the 
reciprocal of the curvature, constant curvatures yield constant radii of curvature, 
and therefore 

R = 1//C(j constant (lO) 

The governing equation for warpage corner deflection, for example, due to the 
curvature along the y-direction is obtained from the geometric relationships depicted 
in figure 2. Thus 
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Q » Sin ^(b/R ) 

y y 


( 12 ) 


6 = R (1 - cos 9 ) 

y y' y 


( 13 ) 


Similarly, for the warpage due to curvature along the x-direction: 


“ ^Acxx 


(ii») 


0 sin”^(a/R ) 


( 15 ) 




(16) 


The comer deflection due to twisting of Boint C, figure 1, is determined 
from recalling that twisting takes place about an axis throu^ the Laminate center: 




( 17 ) 


d » sin ^(a/2R ) 


( 18 ) 
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5 = (1 - cos 0 ) 

^xy xy xy^ 

For twisting on the other edge, a is replaced by b in equation (l8). The total 
corner deflection at point C, figure 1, is obtained by superposition since Unear 
thin-laainate bending theory is xxsed in deriving the curvatures. This procedure 
achieves the objective of obtaining an approximate solution for the corner deflec- 
tions wlthovtt having to solve the coupled fourth-order partial differential equations. 

In addition to the approximations inherent in the Unear thin-laminate 
bending theory, two other approximations were introduced in deriving equations 
(13), (16), and (19). These are: 

1, The projected deformed laminate edge length is approximately equal to the 
undeformed edge length. 

2. Warping due to twisting takes, place about an axis throu^ the laminate 
center. 

Using these approximations, the laminate warpage corner deflections are readily 
derived when the local curvatvures are known. The error introduced by the two 
approxifflations indicated above is not known. 

THEORETICAL WARIA® CORNER DEFLECTIONS 

In this section, theoretically csd-culated data sure presented for warpage comer 
deflections resulting from smaU perturbations of the bending synmietry tislng the 
equations developed in the previous section. The predicted values were obtained 
by perturbing the ply orientation angles in the laminates (0^+30, -30, 0^) and 
(02,+45,-45,Og), using the computer code (ref. U) and the foUowing procedure: 


1 . 


Pertiarb the ply orientation angles in the laminates ( 0^+30, -30, 0^) 
and 

2. Calculate the local curvatures je , >, andjc due to the perturbations 

cxx cyy cxy 

using the comiJuter code in reference 4. 

3. Calculate the corresponding corner deflections using equations (13), 

(l6), and (l9). 

The results obtained by this procedure su’e summarized in table 1, 

•.n this and subsequent tables, the ].suninate configuration is given in the first 

column. The next three colusins are the local curvatures computed using the computer 

code (ref. 4). The corner deflection due to-sc is computed using equation (l3) 

cyy 

and that due using equation (l9). The maximvan warpage comer deflection is 

shown in the last column and is the algebraic sum of the two previous columns. 

The values for x are negligible relative to x and as a result are not included 

cyy 

in the corner deflection calculations. 

A better assessment is obtained by plotting tip deflections versus ply 
orientation perturbations. This is shown in figure 3 for the laminate (Og,+30,-30,0g). 
The following observations can be made: 

1. The warpage corner deflections vary linearly with the perturbation angle. 

2. A perturbation angle of +1.35°, for example, in the +30° pl^ is sufficient 
to produce a warpage tip deflection of .i-,2 inches which is approximately 
equal to four times the laminate thickness which, in this case, is 0.060 
inches . 

A corresponding plot for the {0^+h3,-h3y0^) laminate is shown In figure 
4, A corner deflection of 1.20 inches, for eocample, is obtained at a perturbation 
angle of 11.5 of one of the +45° plies. This exercise shows that it is possible 


to obt»ii* a corner deflection 20 times the laminate thickness with a 11. 5 error 
in only one of the plies. Note the small nonlinearity in the c\irves in figure 4. 

'fhe effects on the corner deflection of replacing one of the -45° plies 
with a +« ply are shown in figure 5. As can he seen in this figure, an angle 9 
of ji'bout l8° will produce a corner deflection of 1.20 inches. The effecwS on 
the corner deflection of the simultaneous perturbation in both the "^45 and —45 
pUes are shown in figure 6. These types of perturbations produce negligible corner 
deflection due to twisting. However, they produce substantial comer deflection 
due to bending. 

The ia^rtant observations firom the previous discussion are as follows: 

1. The warpage corner deflections vary approximately linearly with small 
perturbations in only one of the plies. 

2. Relatively small ply misorientations in only one ply can cause comer 
deflections several times the laminate thickness. 

3 . Conibinations of I'erturbations may be selected to produce specific comer 
deflections . 


SENSITIVITY .'USALYSIS 

In this section, calculated results are presented which can be used to assess 
the significance of the following factors on laminate warpage: 

1. Bending nonsymmetry (two ply equivalent laminate) 

2 . Ply stacking seou nee 

3 . Fiber content nonuniformity (assinning voids in some plies) 


4- Conibinations of items 2 and 3 above 
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The results for the bending nonsymnetry case are summarized In table 2. 

It is seen from the results in this table that only the twisting component is present 
in this case. The graphical representation of the data in table 2 is shown in 
figure 7. Here it is seen that the corner deflection varies nonlinear ly with ply 
angle. It increases rapidly with increasing ply angle in the range 9 10° 

and levels off at 9 >■ 20°. 

Results for some additional ply stacking sequence cases (not examined previously) 
are sumoMrlzed in table 3. As can be seen from the results in this table, ply 
stacking sequences can be selected to yield corner deflections 40 times the laminate 
thickness. Note also that laminates can be made with ply stacking sequences which 
have (l) approodmately the same bending and twisting curvatures and (?) identical 

eurvatiires in the x suid y directions. 

The resvilts from the fiber content nonuniformity sensitivity analysis are 
sumnarized in table h. For this case only, the laminate ( 0^+45, ;U5,Og) was used 
and the fiber content nontmiformlty was introduced via voids. As can be seen from 
the resTilts in table 4, fiber content nonuniformity in the form of about five percent 
voids in only two plies produces corner deflection approximately equal to the 
laminate thickness. Note that fiber content nonxaniformity of the type investigated 
here causes relatively siucU curvature in the x-direction as compared with the other 
two. The graihiical representation of the results in table 4 are shown in figure 8 
when the fiber content nonuniformity is expressed as a void voltmie ratio. As can be 
seen fl-om the curves in this figure, the corner deflections vary nonlinearly with 
fiber content nonunlformi*-y (void volume ratio) and Increases more rapidly as the 
ainotmt of nonuniformity becomes large. Note that at void volume ratios of .10 or 
less, the curves in figin*© 8 are approxiinately Linear* 
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Resiilts of sensitivity analysis for the comhined case of small ply misorienta- 
tions and fiber content nonuni formity are summarised in table 5. As can be seen 
fraa the data in this table, the combinations selected, which on the surface appear 
to be minor variations from case to case, produce substantially different corner 
deflections. Comparing corner deflections from the individual cases in tables 1 
and k with the corresponding values in table 5 , it is seen that some combinations 
are compensatory. The tabvilated curvatures in tables 1-5 can be used to predict 
laminate varpage corner deflections in laminates of similar cmifiguration but 
different planform dimensions. For these cases, the edge dimensions are the only 
variables that change in equations (u)-(l9)> 

It is noted here that the same procedure used for the corner deflections can 
be used to calculate varpage deflection at any point (x, 39 fig. l) in the laminate. 
For these calculations, the sign of the subtended angle is detennlned from the 
foUowing equations: 

= sin“^(y/Ey.) (20) 

= sin“^(x/E^) (2l) 



sin‘^(x/2R ) 


(22) 




(23) 
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where the physical problem will dictate whether equation (22) or (23), or some 
ccanbination thereof, should be used. 

The la5)ortant points from the previous discussion are: 

1. Relatively small perturbations in various ply stacking sequences and fiber 
content nonundformitles produce considerable deflections, 
t 2. At small values of the perturbation, the corner deflection varies approxl* 

j mately linearly with the perturbation. 

! 

3, Conblned perturbations can produce compensatory effects on the corner 
deflection compared to that produced by the Individual cases. 

APPLICATION OF EQUATIONS FOR CALCULATING WARBAGE 

The apprarimate equations for warpage corner deflections derived previously 
can be \ised to sxiggest possible combinations of nonsymnetrles that cotild explain 
the warpage measured in two actual laminates. 

Description of Laminates 

Laminates were fabricated to have symmetry with respect to bending and 

i 

I siqiplled to NASA-Lewls Research Center by an outside vendor. These laminates were 
I 12 inch square plates, eight plies thick (O.O6O inches) of ply 

I stacking sequence, and were made from MOD-l/ERLA U6l7 with about 50 i>ercent fiber 

j 

I volume ratio and cured at about 370°F. Two of these lamlisates, (0g,+30,-30,0g) 

5 

I and (Og+U5,-45,0-)^ exhibited warpage with measurable corner deflections as depicted 

i in figure 1. After curing, the lamixiates were reheated in an unrestrained condition 

f 

I to about 300 F and then cooled to room temperature. 
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Corner Deflection Measurements 

The varpage corner deflections of the tvo laminates (O^U5,-U5,Og) and 
(0_,+30,-30>'*‘*»^5) were measiired at points C and D as depicted in figure 1. 

The measured values are summarised in table 6, Note that the last tvo columns in 
table 6 separate the corner deflections into tvo components: bending and tvisting. 

The bending component eqvials the algebraic average of the tvo corner deflections 
and the tvisting component eqvials their algebraic difference. . As is observed from 
the values shovn in table 6, the maxlnnnn corner deflections can be substantial: 

1.20 inches for the (02,+U5,;U5,Og) laminate and 0.22 inches for the (0^30, +30,0g) 
laminate. Cospared to 0.06 iziches for the lazoinate thickness, these deflections are 
20 times the laminate thickness for the laminate and about four 

times for the ( 0^+30, -30, 0^) laminate. 

Scsne Bossible Combinations of Nonsynsnetrles 

The foUovlng procedure vas used to identify a fev of a large number of possible 
combinations of nonsymnetries that could have contributed to the warpaje corner 
deflections that vere measured in the actiial laminates. 

Fat the (0^+30,-30,02) laminate, the data in table 1 were used as a guide 
hjecause the 12- inch edge of the actual laminate is only 20 percent longer than the 
10-inch edge which vas used to generate the ctirves in figure 3* Comparing corres- 
|)ondlng values from table 6 and table 1, it is seen that the measured value of the 
m^T<nniw corner deflect!' b^'tveen those calculated in the lamdnates 

(Og,31,-30,;^30, Og) and vO^, ' 3^,-30;30,02). Using a=b.»12 in equations (ll) to (l3) 
and (17) to (19) and plotting the data as shown in figure 3 yields a perturbation 
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nonsymnetry of 0.9°. The components of the maximian warpage deflection produced by 
this pertvjrbation are: 

0 .I 7 U inches due to bending (.17 inches measured) 

0.046 inches due to twisting (.05 inches measured) 

It is of interest that the measured values are quite similar to those calcu- 
lated for the arbitrarily assumed perturbation. 

If the ply orientations in the laminate were (0^45,;39.5,0g), 

the calculated amximum warpage corner deflection would be 1.20 inches which is 
equal to that measured. 

The corresponding components of the deflection are: 

1,175 inches due to bending (l.OO inches measured) 

0.025 Inches due to twisting (0.20 inches measured) 

Agaiiv it is of interest to note that the measured values are quite similar 
to those calculated for the arbitrarily assxamed perturbation and that the perturba- 
tions required to produce the warpage corner deflections are relatively small. 

GENERAL RMAEKS 

The equations derived herein together with the results presented should 
be of considerable inractical value to both designers and fabricators of fiber composites. 
For example, the fabricator can use the equations to obtain an a priori assessment 
of the warpage resulting from fskctors which may be difficult to control accurately 
d\irlng laminate fabrication. The designer may use the results to specify tolerances 
which will minimiee warpage-producing nonsymmetries; or he may use the results to 
design the laminate with warpage-compensating nonsymmetries. 


The results of this investigation also illustrate that suitable combinations 
of warpage-producing factors are compensatory and can be used to alleviate some 
fabrication problems. It is also noted that because a laminate is flat, or meets 
some flatness specification tolerances, it may not be free of warpage-producing 
nonsymmetries. It is possible that warpage will result if the laminate is reheated 
and cooled again. This type of warpage may cause difficulties in designs where the 
laminate is part of a stinicture restrained along its edges and is subjected to 
thermal cyclic environment. In this case, alternating stresses due to warpage would 
not have been taken into account in designing the laminate. 

Althou^ the data presented were not obtained using hybrid composites, such 
composites will exhibit similar warpage in the presence of bending nonsynmetries 
and residual stresses. The approximate equations derlged herein are applicable to 
these composites as well. 


SIMiARY OF RESULTS 

The major results of this investigation are s\mBnarized below: 

1. A convenient set of equations has been derived which can be used to 
approximate laminate warpage corner deflections resulting from small 
errors during fabrication. This approach does not require solutions of 
coupled fourth-order partial differential equations which would normally 
be the case. 

2. Using the derived equations, it was found that an error of only one 
degree in one ply orientation was sufficient to produce a warpage corner 
deflection of a magnitude twice the laminate thickness. This ply mis- 
orientation, among a large number of other possible ones, could have 
caused the warpage corner deflection that was measured in a 10 inch x 
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10 inch X .06 Inch (Q^;^30,-30»92^ laminate made from MOD-I/ERLA-U617. 

3. Using the derived equations in a sensitivity anaj^sis of the various 
factors contributing to laminate warpage, it was found that a 3^^ 
fiber migration or a five percent void volime ratio in some plies is 
sufficient to produce laminate warpage corner deflections several times 
the ply thickness. It was also found that errors in ply stacking sequences 
can produce laminate warpage corner deflections as much as Uo times the 
laminate thickness . 

U. Combinations of errors or tolerances can produce compensating effects 
on the warpage and initially flat laminates may warp as a result of 
unrestrained reheating and cooling. 

5. The warpage corner deflection varies linearly with errors up to 10° 

in ply misorientations and vtp to 10 percent in fiber content nonuniformity. 


I 
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SYMBOLS 

Array of composite axial stiffnesses referred to composite structural 
axes 

Laminate edge dimension along x 
Laminate edge dimension along y 

Array of composite coupling stiffnesses referred to composite 
structural axes 

Array of eoaqposite bending (flexural) stiffnesses referred to composite 
structural axes 

Array of strain-stress relations (elastic constants for the ply) 

Ply longitudinal modulus 
Ply transverse modulus 
Ply shear modulus 

Fiber and void volume ratios, respectively 

Vector of unbalanced thermal moments referred to composite structural 
axes 

Vector of unbalanced thermal forces referred to composite structural 
axes 

Radius of curvature (subscripts denote direction) 

Array of transformation coefficients for the i*'^ ply 
Temperature difference for the 1^^ ply 
Structural axes 

Location of the i^^ ply relative to the reference plane 
Material axes 



®^i22 


^cox 

V|12 

^t21 


Vector of thermal coefficients of expansion of the i^'' ply 
Ply longitudinal thermal coefficient of expansion 
Ply transverse thermal coefficient of expansion 

Corner deflection (subscript denotes compo /int due ^o corresponding 
CTorvature ) 

Vector of composite strains referred to coamposite structural axes 
at the reference plane 

Ply angle meas\ired from the ccssposite structural axes to the ply 
material axes 

Vector of composite local curvatures referred ro cosiposite structxiral 
axes 

Fly major Poisson's ratio 
Ply minor Poisson's ratio 


cxx 

Local 

ctirvature 

•^cyy 

Local 

curvatture 

*^cxy 

Local 

curvature 
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TABLE 1 - THEORETICALLY CALCULATED EFFECTS OF SMALL PLY MISORI- 
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^Here and In subse<|uent tables unless otherwise noted: 

1. The composite material is MODMOR-l/ERLA 4617 at . 5 fiber volume ratio 
and cured at 870° F yielding a AT of -300° F. 

2. The laminate dimension are a = b = 10 in, : h = . 06 in. (fig. 1). 








TABLE 2. - THEORETICALLY CALCULATED EFFECTS OF BENDING 
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TABLE 6. - MEASURED PLATE WARPAGE CORNER DE- 


FLECTIONS (REF. FIGURE 1, a=b=12 IN. , THICKNESS = 
.06 IN. GRAPHITE FIBER MOD- I/EPOXY ERLA 4617) 


Laminate 

Corner deflections, in. 


Maximum 

Minimum 

Due to 


p<^ir.t C 

point D 

bending^ 

twisting 

[Oj.rfS, =t45,02] 

1.20 

.80 

1.00 

±.L0 

[02,±30,t30,0^] 

.22 

.12 

.17 

t.05 


^Bending component = 1/2 (deflection at C + deflection at D) 
“Twisting component =1/2 (deflection at C - deflection at D) 
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f i^urf 1. - SchcmaUc of (telormed nonsymimlrlc anglepHfd lamirolt with residual 
stre>se$. 





Figure 2. - Set/ T«tk: depicting geofnetric relationships between 
corner deflection vnel ''ilge (finienston and raefius of curvature. 





COMBINED, 



Figure 5. - Predicted maximum warpage corner 
deflection (point C, fig. 1, a - b - 10 in. ) due 
to residual stress in a [(k. *45, -45+A0, +45,021 
laminate from MOD- I/ERiA *4617 con^ite at 
0. 5 fiber volume ratio and AT • -300° F. 



Figure 6. - Predicted maximum warpage corner 
deflection (point C, fig. 1, a > b * 10 in. l due 
to residual stress in a lOo. ±45.^5FAO,0^ 
laminate from MOD-I/ERLA 4617 composite at 
0.5 fiber volume ratio and AT • -30(r F, 
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TWIST m BENDING IN 
THIS CONFIGURATION) 





.10 ,20 .30 
VOID VOLUME RATIO, 


Figure 8. - Predicted maximum warpage corner deflection (point C. 
fig. 1. a - b - 10 in. ) due to residual stress in a 
laminate with void volume ratio ky 14(0,0). 4(i(y)l made from 
MOO-I/ERLA-4617 composite at 0.5 fiber volume ratio and AT • 
-300® F, 
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